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Abstract 

The overview is given of the results obtained recently in the course of renormalization-group 
(RG) study of two-dimensional (2D) models. RG functions of the two-dimensional n-vector Xcj) 4 
Euclidean field theory are written down up to the five-loop terms and perturbative series are 
resummed by the Pade-Borel-Leroy techniques. An account for the five-loop term is shown to shift 
the Wilson fixed point only briefly, leaving it outside the segment formed by the results of the 
lattice calculations. This is argued to reflect the influence of the non-analytical contribution to the 
/3-function. The evaluation of the critical exponents for n = 1, n = and n = — 1 in the five-loop 
approximation and comparison of the results with known exact values confirm the conclusion that 
non-analytical contributions are visible in two dimensions. The estimates obtained on the base of 
pseudo-e expansions originating from the 5-loop 2D RG series are also discussed. 
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The field-theoretical renormalizat ion-group (RG) approach proved to be a powerful tool 
for calculating the critical exponents and other universal quantities of the basic three- 
dimensional (3D) models of phase transitions. Today, many-loop RG expansions for f3- 
functions (six-loop), critical exponents (seven-loop), higher-order couplings (four-loop), etc. 
of the 3D 0(n)-symmetric, cubic, and some other models are known resulting in high- 
precision numerical estimates for experimentally accessible quantities [1-7]. The main aim 
of this paper is to demonstrate how effective (or ineffective) is the field-theoretical RG ma- 
chinery in two dimensions where i) the RG series are stronger divergent and ii) singular 
(non- analytic) contributions to RG functions are expected to be larger than for 3D systems. 

The Hamiltonian of the model describing the critical behavior of various 2D systems 
reads: 



H= d 2 x 



^(m^ + (V^) 2 ) + A(^ )2 



(1) 



where (p a is a real n- vector field, m 2 , is proportional to T — TjP\ Tj°) being the mean- field 
transition temperature. 

The /9-function and the critical exponents for the model (1) are calculated within the 
massive theory, with the Green function, the four-point vertex and the <p 2 insertion being 
normalized in a conventional way: 

dGR 1 (p,m,g4) 



(■!,■■ ((I II). Hi) = V) 



2 



= 1, (2) 



dp 

r fl (0, 0, 0, m, g) = m 2 g 4 , rf(0, 0, m, g 4 ) = 1. 

Since the four-loop RG expansions at n — 1 have been obtained many years ago [1], we 
are in a position to find corresponding series for arbitrary n and to calculate the five-loop 
terms. The results of our calculations are as follows [8]: 
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(n + 8) 5 

the refined expression for the five-loop contribution to 7 _1 is taken from [9]. Instead of the 
renormalized coupling constant g^, a rescaled coupling 

9 = ~2A^ 9A ' (6) 

is used as an argument in above RG series. This variable is more convenient since it does 
not go to zero under n — > oo but approaches the finite value equal to unity. 

To evaluate the Wilson fixed point location g* and numerical values of the critical expo- 
nents, the resummation procedure based on the Borel-Leroy transformation 

oo °9 oo 

f(x) = £ Cl x l = J e-H b F(xt)dt, F(y) = ^ 7-7W » ( 7 ) 



i=0 







is used. The analytical extension of the Borel transforms is performed by exploiting relevant 
Pade approximants [L/M]. In particular, four subsequent diagonal and near-diagonal ap- 
proximants [1/1], [2/1], [2/2], and [3/2] turn out to lead to numerical estimates for g* which 
rapidly converge, via damped oscillations, to the asymptotic values; this is cleary seen from 
Table 1. These asymptotic values, i. e. the final five-loop RG estimates for g* are presented 
in Table 2 for < n < 8 (to avoid confusions, let us note that models with n > 2 possessing 
no ordered phase are studied here only as polygons for testing the numerical power of the 
perturbative RG technique). As Table 2 demonstrates, the numbers obtained differ appre- 
ciably from numerical estimates for g* given by the lattice and Monte Carlo calculations 
[10-16]; such estimates are usually extracted from the data obtained for the linear (x) and 
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non-linear (x*) susceptibilities related to each another via g±: 



d 3 M 2 _ 2 

X4 = "oTjT =~xm 04, 



(8) 



H=0 



Since the convergence of the numerical estimates for g* given by the resummed RG series is 
oscillatory, an account for higher-order (six-loop, seven-loop, etc.) terms in the expansion 
(3) will not avoid this discrepancy [9]. That is why we believe that it reflects the influence 
of the singular (non-analytical) contribution to the /^-function. 

The critical exponents for the Ising model (n — 1) and for those with n = and n = — 1 
are estimated by the Pade-Borel summation of the five-loop expansions (4), (5) for 7 _1 and 
77. Both the five-loop RG (Table 1) and the lattice (Table 2) estimates for g* are used in 
the course of the critical exponent evaluation. To get an idea about an accuracy of the 
numerical results obtained the exponents are estimated using different Pade approximants, 
under various values of the shift parameter b, etc. In particular, the exponent 77 is estimates in 
two principally different ways: by direct summation of the series (5) and via the resummation 
of RG expansions for exponents 



which possess a regular structure favouring the rapid convergence of the iteration procedure. 
The typical error bar thus found is about 0.05. 

The results obtained are collected in Table 3. As is seen, for small exponent rj and in some 
other cases the differences between the five-loop RG estimates and known exact values of the 
critical exponents exceed the error bar mentioned. Moreover, in the five-loop approximation 
the correction-to-scaling exponent uj of the 2D Ising model is found to be close to the value 
4/3 predicted by the conformal theory [17] and to the estimate 1.35 ± 0.25 extracted from 
the high-temperature expansions [18] but differs markedly from the exact value uj = 1 [19] 
and contradicts to the conjecture uj = 2 [20]. This may be considered as an argument in 
favour of the conclusion that non-analytical contributions are visible in two dimensions. 

The field theory enables us also to find the higher-order, sextic coupling constant entering 
the free energy expansion in powers of the magnetization M. For 2D Ising model 



where m is a renormalized mass, Z being a field renormalization constant. In the critical 
region, where fluctuations are so strong that they completely screen out the initial (bare) 
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interaction, the behaviour of the system becomes universal and dimensionless effective cou- 
plings gik approach their asymptotic limits g\ k . 

In order to estimate g\ we calculate RG expansion for g 6 and then apply Pade-Borel- 
Leroy resummation technique to get proper numerical results. As is well known, accurate 
enough numerical estimates may be extracted only from sufficiently long RG series. We 
have obtained the expression for g 6 in the four-loop approximation [21] which turned out to 
provide fair numerical estimates for the quantity of interest. 

The method of calculating the RG series we used in [21] is straightforward. Since in two 
dimensions higher-order bare couplings are irrelevant in RG sense, renormalized perturbative 
series for g e can be obtained from conventional Feynman graph expansion of this quantity in 
terms of the only bare coupling constant - quartic coupling A. In its turn, A may be expressed 
perturbatively as a function of renormalized dimensionless quartic coupling constant g±. 
Substituting corresponding power series for A into original expansion we can obtain the RG 
series for g 6 . As was shown [22,23], the one-, two-, three- and four-loop contributions are 
formed by 1, 3, 16, and 94 one-particle irreducible Feynman graphs, respectively. Their 
calculation along with the renormalization procedure just described gives: 

g 6 = — g \ (\ - 3.2234882 g 4 + 14.957539 g\ - 85.7810 gf). (11) 

This series may be used for estimation of the universal number g$. 

With the four-loop expansion in hand, we can construct three different Pade approx- 
imants: [2/1], [1/2], and [0/3]. To obtain proper approximation schemes, however, only 
diagonal [L/L] and near-diagonal Pade approximants should be employed. That's why fur- 
ther we limit ourselves with approximants [2/1] and [1/2]. Moreover, the diagonal Pade 
approximant [1/1] is also dealt with although this corresponds, in fact, to the usage of the 
lower-order, three-loop RG approximation. 

Since the Taylor expansion for the free energy contains as coefficients the ratios R 2 k = 
92k/ dt~ l rather than the renormalized coupling constants themselves: 

F(z) - F(0) = ^(j + z 4 + ^ + Rsz 8 + .-), z 2 = 9 -^f, (12) 

we work with the RG series for i? 6 . It is resummed in three different ways based on the 
Borel-Leroy transformation and the Pade approximants just mentioned. The Borel-Leroy 
integral is evaluated as a function of the parameter b under g 4 = g\. For the fixed point 
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coordinate the value g\ = 0.6124 [24] is adopted which is believed to be the most accurate 
estimate for g\ available nowadays. The optimal value of b providing the fastest convergence 
of the iteration scheme is then determined. It is deduced from the condition that the Pade 
approximants employed should give, for b = b opt , the values of R* & which are as close as 
possible to each other. Finally, the average over three estimates for is found and claimed 
to be a numerical value of this universal ratio. 

The results of our calculations are presented in Table 4. As one can see, for b = 1.24 all 
three working approximants lead to practically identical values of i?g. Hence, we conclude 
that for 2D Ising model at criticality 

R* 6 = 2.94, gl = 1.10. (13) 

How close to their exact counterparts may these numbers be? To clear up this point let 
us discuss the sensitivity of numerical estimates given by RG expansion (11) to the type of 
resummation. The content of Table 4 implies that, among others, the results given by Pade 
approximant [2/1] turn out to be most strongly dependent on the parameter b. This situation 
resembles that for 3D 0(n) -symmetric model where Pade approximants of [L — 1/1] type 
for /5-function and critical exponents lead to numerical estimates demonstrating appreciable 
variation with b while for diagonal and near-diagonal approximants the dependence of the 
results on the shift parameter is practically absent [1,3,25]. In our case, Pade approximants 
[1/1] and [1/2] may be referred to as generating such "stable" approximations for g$. Since 
for b varying from to 15 (i.e., for any reasonable b) the magnitude of gl averaged over these 
two approximations remains within the segment (1.044, 1.142) it is hardly believed that the 
values (13) can differ from the exact ones by more than 5%. 

Another way to judge how accurate our numerical results are is based on the comparison 
of the values of gl given by four subsequent RG approximations available. While within the 
one-loop order we get gl = 2.633 which is obviously very bad estimate, taking into account of 
higher-order RG contributions to g 6 improves the situation markedly. Indeed, two-, three-, 
and four-loop RG series when resummed by means of the Pade-Borel technique with use of 
"most stable" approximants [0/1], [1/1], and [1/2] yield for gl the values 0.981, 1.129, and 
1.051, respectively. Since this set of numbers demonstrates an oscillatory convergence one 
may expect that the exact value of renormalized sextic coupling constant lies between the 
higher-order - three-loop and four-loop - estimates. It means that the deviation of numbers 
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(13) from the exact values would not exceed 0.05. 

It is instructive to compare our estimates with those obtained by other methods. S.-Y. 
Zinn, S.-N. Lai, and M. E. Fisher analyzing high temperature series for various 2D Ising 
lattices found that Rq = 2.943 ±0.007 [26]; almost identical value was obtained in [27]. Our 
result for Rq is seen to be in a brilliant agreement with this number. Of course, practical 
coincidence of the lattice and four-loop RG estimates is occasional and can not be considered 
as a manifestation of extremely high accuracy of the methods discussed. The closeness of 
these estimates to each another, however, unambiguously demonstrates high power of both 
approaches. Moreover, such a closeness shed a light on the role of a singular contribution to 
<?6 which can not be found perturbatively: this contribution is seen to be numerically small. 

It is interesting also to address the results given by another field-theoretical approach - 
the e expansion. For the Ising systems three terms in the e expansion for i? 6 are known [28]: 



Let us apply a simple Pade-Borel procedure to this series as a whole and to the series in 
brackets and then put e = 2. We find R* & = 3.19 and R* & = 3.12 respectively, i.e. the numbers 
which differ from our estimate by less that 9%. Keeping in mind lack of a small parameter 
these values of R§ may be referred to as consistent. Proper account for higher-order terms 
in the e expansion for R e should make corresponding numerical estimates closer to those 
extracted from 2D RG and high-temperature series. Very good agreement between the first 
number (13) and the estimate R* Q = 2.95 ±0.03 [27] obtained by matching of the e expansion 
with the exact results known for D = 1 and D = may be considered as an argument in 
favor of this belief. One can find more details in recent comprehensive review [29]. 

Along with the RG calculations at physical dimension and the e expansion, some other 
field-theoretical approach may be employed to estimate the critical parameters of 2D Ising 
model. We mean the method of the pseudo-e expansion (see Ref. 19 in [2]). Pseudo-e 
expansions for the Wilson fixed point coordinate g* and critical exponents can be easily 
derived from the RG series (3)-(5) using standart technique. They are as follows [30]: 




(14) 



g * = T + 0.716173621r 2 + 0.095042867r 3 + 0.086080396r 4 - 0.204139r 5 , 



(15) 
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7] = 0.033966147r 2 + 0.046628762r 3 + 0.030925471r 4 + 0.0256843r 5 . (17) 

Note that the higher-order terms in series (15) and (16) have small numerical coefficients 
and irregular signs. Smallness of these coefficients enables one to obtain accurate enough 
estimates for g* and critical exponent 7 without addressing the Borel- transformation-based 
resummation methods. 

To demonstrate this, conventional Pade triangles originating from (15) and (16) under r = 
1 are presented here (Tables 5 and 6). Since diagonal and near-diagonal Pade approximants 
are known to exhibite the best approximating properties, the numbers 1.751 and 1.837 from 
Table 5 should be referred to as most reliable estimates for g*. Averaging over them, we 
obtain g* = 1.794 which differs from the exact value g* = 1.75436 [24] by 2%. As seen from 
Table 5, it is the five- loop approximation that provides so good numerical result; almost 
all lower-order approximations suffers from dangerous poles resulting in strongly scattered 
estimates. The same is true for the susceptibility exponent. Indeed, the numbers given 
by the main working approximants [2/3] and [3/2], as well as by approximant [4/1], are 
almost coincide with each other and are close to the exact value 7 = 1.75. In contrast, 
approximants [2/2] and [1/3], corresponding to the four-loop order, have dangerous poles 
which considerably affect the results. 

Unfortunately, the pseudo-e expansion technique turns out to be much less powerful when 
applyed to estimate "small" critical exponent 77. Both the direct summation of the expansion 
(17) and Pade resummation of the series for "big" exponents 7 and v lead to the numbers 
differing by 0.1 and even more from the exact value rj = 0.25 [30]. To the contrary, the 
pseudo-e expansion for the ratio R 6 = g§jg\ 

R 6 = 4r(l - 0.409036r + 0.305883r 2 - 0.437676r 3 ) (18) 

demonstrates good Pade summability. It is clearly seen from Table 7 [30]. Averaging over 
two working approximants [2/2] and [3/1] gives the number R Q = 2.90 which is close to 
earlier estimates R 6 = 2.94 [21], R 6 = 2.95 [27], R 6 = 2.943 [26], and to high-precision 
values R 6 = 2.94294 [24], R 6 = 2.94238 [29,31]. Usage of more advanced, Pade-Borel 
resummation technique shifts the pseudo-e expansion estimate to Re = 2.94 [30] making it 
practically equal to just mentioned numbers. 
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The area where 2D A0 4 field theory can be successfully applyed is not limited by Ising- 
like and 0(n) -symmetric systems. The RG analysis of 2D cubic, MN, chiral, and weakly 
disordered models proofs to be rather effective provided the higher-order - four- and five-loop 
- approximations are used [9,32-34]. In particular, many-loop RG calculations reproduce 
with high accuracy the exact results known for 2D anisotropic systems with n-vector order 
parameters. Detailed description of the situation may be found in [9,33,34]. 

I thank P. Calabrese, D. V. Pakhnin, P. Parruccini, and E. V. Orlov for fruitful collabo- 
ration. I am also grateful to B. N. Shalaev for numerous valuable discussions of the critical 
thermodynamics of 2D systems. This work was supported by the Russian Foundation for 
Basic Research under Grant No. 04-02-16189. 
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TABLE I: The Wilson fixed point coordinate for models with n = 1, n = and n = — 1 in four 
subsequent RG approximations and the final five- loop estimates for g*(n). 



n 


[1/1] 


[2/1] 


[2/2] 


[3/2] 


g*, 5-loop 


1 


2.4246 


1.7508 


1.8453 


1.8286 


1.837 ± 0.03 





2.5431 


1.7587 


1.8743 


1.8402 


1.86 ± 0.04 


-1 


2.6178 


1.7353 


1.8758 


1.8278 


1.85 ± 0.05 



TABLE II: The Wilson fixed point coordinate g* and critical exponent u> for < n < 8 obtained in 
the five-loop RG approximation. The values of g* extracted from high-temperature (HT) [11,13] 
and strong coupling (SC) [12] expansions, found by Monte Carlo simulations (MC) [14,15], obtained 
by the constrained resummation of the e-expansion for g* (e-exp.) [13], and given by corresponding 
1/n-expansion (1/n-exp.) [13] are also presented for comparison. 



n 





1 


2 


3 


4 


8 


9* 


RG, 5-loop 


1.86(4) 


1.837(30) 


1.80(3) 


1.75(2) 


1.70(2) 
(6=1) 


1.52(1) 
(6 = 1) 


HT 


1.679(3) 


1.754(1) 


1.81(1) 


1.724(9) 


1.655(16) 




MC 




1.71(12) 


1.76(3) 


1.73(3) 






SC 


1.673(8) 


1.746(8) 


1.81(2) 


1.73(4) 






e-exp. 


1.69(7) 


1.75(5) 


1.79(3) 


1.72(2) 


1.64(2) 


1.45(2) 


1/n-exp. 








1.758 


1.698 


1.479 


to 


RG, 5-loop 


1.31(3) 


1.31(3) 


1.32(3) 


1.33(2) 


1.37(3) 


1.50(2) 
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TABLE III: Critical exponents for n = 1, n = 0, and n = — 1 obtained via the Pade-Borel summa- 
tion of the five-loop RG expansions for 7 -1 and 77. The known exact values of these exponents are 
presented for comparison. 



n 




* 

9 


7 


V 


V 


a 


1 


RG 


1.837 


1.79 


0.146 


0.96 


0.07 






1.754 (HT) 


1.74 


0.131 


0.93 


0.14 




exact 




7/4 


1/4 


1 











(1.75) 


(0.25) 









RG 


1.86 


1.45 


0.128 


0.77 


0.45 






1.679 (HT) 


1.40 


0.101 


0.74 


0.52 




exact 




43/32 


5/24 


3/4 


1/2 








(1.34375) 


(0.20833) 


(0.75) 


(0.5) 


-1 


RG 


1.85 


1.18 


0.082 


0.62 


0.76 






1.473 (SC) 


1.15 


0.049 


0.59 


0.82 




exact 




37/32 


3/20 


5/8 


3/4 








(1.15625) 


(0.15) 


(0.625) 


(0.75) 



TABLE IV: The values of Rq obtained by means of the Pade-Borel-Leroy technique for various b 
within three-loop (approximant [1/1]) and four-loop (approximants [1/2] and [2/1]) RG approxima- 
tions. The estimate for b = 1 in the middle line is absent because corresponding Pade approximant 
turnes out to be spoilt by a positive axis pole. 

b 1 1.24 2 3 4 5 7 10 15 
[1/1] 2.741 2.908 2.937 3.009 3.077 3.125 3.161 3.212 3.258 3.301 
[1/2] 2.827 - 2.936 2.877 2.853 2.838 2.828 2.814 2.800 2.787 
[2/1] 3.270 2.988 2.936 2.800 2.667 2.568 2.491 2.380 2.273 2.171 
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TABLE V: The Wilson fixed point location g* , extracted from pseudo-e expansion (15) by means of 
constructing Pade approximants [L/M]. Coordinates of " dangerous" poles of Pade approximants, 
i. e. those lying on the real positive semiaxis are indicated as subscripts. 



M L 


1 


2 


3 


4 


5 





1.000 


1.716 


1.811 


1.897 


1.693 


1 


3.523i.4 


I.8267.5 


2.724i.i 


1.837 




2 


1.425 


I.9I83.0 


1.850 6 .i 






3 


2.6OI1.4 


1.751 








4 


1.194 











TABLE VI: Numerical values of the critical exponent 7 obtained by Pade summation of series (16) 
for 7 -1 . 



M L 





1 


2 


3 


4 


5 





1.000 


1.500 


1.808 


1.730 


1.859 


1.660 


1 


1.333 


2.024 2 . 9 


1.744 


1.778 


1.777 




2 


1.558 


1.702 


I.8OO5.2 


1.777 






3 


1.646 


6.871i.i 


1.772 








4 


1.732 


1.718 










5 


1.7146.1 













TABLE VII: Pade triangle for the universal ratio Rq given by pseudo-e expansion (18). 



M L 


1 


2 


3 


4 





4.000 


2.364 


3.587 


1.837 


1 


2.839 


3.064 


2.867 




2 


3.148 4 .5 


2.940 






3 


2.621 
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